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Using the tensor network approach, we investigate the monomer-dimer models on a checkerboard lattice, in 
which there are interactions (with strength v) between the parallel dimers on one-half of the plaquettes. For 
the fully-packed interacting dimer model, we observe a Kosterlitz-Thouless (KT) transition between the low- 
temperature symmetry breaking and the high-temperature critical phases; for the doped monomer-dimer case 
with finite chemical potential p, we also find an order-disorder phase transition which is of second-order, instead. 

We use the boundary matrix product state approach to detect the KT and second-order phase transitions, and 
obtain the phase diagrams v — T and p — T. Moreover, for the non-interacting monomer-dimer model (setting 
p = y = 0), we get an extraordinarily accurate determination of the free energy per site (negative of the 
monomer-dimer constant h 2 } as / = -0.662 798 972 833 746 with the dimer density n = 0.638 123 109 228 547, 
both of 15 correct digits. 
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I. INTRODUCTION 

Classical monomer-dimer model in two dimensions (2D) 
is one of the intriguing models in statistical mechanics. The 
problem has a venerable history,it was firstly introduced 
in the context of the absorption of molecules on the surface.' 
When a rigid molecule occupies two nearest neighbor (NN) 
sites on the (square) lattice, it can be regarded as a dimer link¬ 
ing the two NN sites, while an empty site means the presence 
of a monomer. The monomer-dimer model can be related to 
the Ising^ and the height models,^ playing an important role in 
the statistical physics. A special case of the monomer-dimer 
model, namely the fully-packed dimer model, can be ana¬ 
lytically solved,^’"'’® while the general monomer-dimer case 
is not. Numerating all the possible configurations and cal¬ 
culating the properties of the monomer-dimer model is an 
NP-complete problem and thus unfortunately “intractable” in 
computations.Numerically, one has to adopt some ap¬ 
proximate methods, like the Monte Carlo samplings,to 
study the monomer-dimer models. 

The fully-packed dimer models exhibit different proper¬ 
ties on bipartite and non-bipartite 2D planar lattices. The 
former supports a critical phase with the algebraic decay¬ 
ing dimer-dimer correlations,^’^ while the latter (say, the tri¬ 
angular and kagome lattices) have exponential dimer-dimer 
correlations.Exept the 2D lattices, people have also 
investigated the hard-core dimer models on various 3D 
lattices.'^ Extended critical phases are found on the bipartite 
cubic lattice, while no critical phases are found on the non- 
bipartite 3D lattices.'^ Moreover, the classical dimer models 
can be “upgraded” to the so-called quantum dimer models, 
by promoting the classical dimer configurations to quantum 
state bases. The quantum dimer model was introduced by 
Rokhsar and Kivelson,'* where the singlet formed by two ad¬ 
jacent spins plays the role of a dimer. The quantum dimer 
model is one of the typical systems which exhibit nontrivial 
topology and fractional excitations. 


Recently, E. Alet et al. introduce the interacting dimer mod¬ 
els on a square lattice, where two parallel dimers on the same 
plaquette are coupled (attractive). They studied this interact¬ 
ing dimer model by Monte Carlo simulations, and found a 
dimer order-disorder phase transition, of Kosterlitz-Thouless 
(KT) type,^^ at a certain temperature.'^’'^ Interestingly, the in¬ 
troduction of interactions between dimers is not only of theo¬ 
retical interest in the model study, but also acquire experimen¬ 
tal realizations recently. An adsorption experiment of certain 
rodlike organic molecules on the graphite was reported^^ and 
it is found to be relevant to the fully-packed dimer model on a 
hexagonal lattice, with couplings between neighboring paral¬ 
lel dimers. 


(a) (b) 



FIG. 1: (Color online) (a) The illustration of the monomer-dimer 
model on a checkerboard lattice. Interactions between the parallel 
dimers are introduced on the green (gray) plaquettes. The underlying 
tensor network is also shown, where the plaquette tensors (T) only 
cover one-half of the plaquettes. The square-lattice tensor network 
can be divided into two sublattices, namely A and B, denoted by 
dark and light green (gray) colors, respectively. Besides the plaquette 
tensor T, there is a local matrix P on each vertex, (b) The plaquette 
tensor T is shown, j,- (i 6 j 1,2,3,4j) are the four indices of tensor T. 
(c) The vertex matrix P can be absorbed into a plaquette tensor T by 
contracting the sharing indices (^ 2 ). 
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In this work, we introduce and study the interacting 
monomer-dimer model on a black-white checkerboard lattice. 
We find and employ the compact tensor network (TN) rep¬ 
resentation of the (grand) partition function to investigate the 
model. The problem of calculating the free energy is thus 
transformed into a problem of how to accurately contract the 
TN. In practice, we contract it row by row with the infi¬ 
nite time evolving block decimation (iTEBD) technique for 
the matrix product state (MPS)^^ and calculate the thermo¬ 
dynamic properties of the monomer-dimer model with high 
precision. Through numerical simulations, we show that the 
fully-packed dimer model has a low temperature dimer or¬ 
dered phase and a high temperature critical phase, with a KT 
transition separating these two phases. We observe no singu¬ 
larity in the energy and its derivative (specific heat) curves, 
however, we detect the KT transition by calculating the order 
parameters and correlation length. On the other hand, when 
monomer doping is introduced in the model (with the chem¬ 
ical potential /i), there is also an order-disorder phase transi¬ 
tion at certain temperature, which is instead found to be of 
second order. Remarkably, in addition to the regular thermo¬ 
dynamic quantities including the energy derivatives and the 
order parameter, etc, we also detect the (second-order) phase 
transition by checking the “entanglement” properties of the 
boundary MPS. We even extract the corresponding conformal 
central charge of the high-T critical phase of the fully-packed 
dimer model, by studying the block entanglement entropies 
of the boundary MPS. At last, collecting the phase transition 
points, we show and discuss the jU-T and v-T phase diagrams 
of the interacting monomer-dimer model. 

The rest of this paper is organized as follows. In Sec. II, we 
introduce the TN representation of the partition function and 
the method for accurate evaluation of the thermodynamics. In 
Sec. Ill, we show the main numerical results on phase transi¬ 
tions, and the phase diagrams of the monomer-dimer model. 
The last section (Sec. IV) is devoted to a conclusion of the 
paper. 


II. MODEL AND METHODS 


Setting ju —» oo, we recover the fully-packed dimer model (no 
monomer doping). 

The partition function of the monomer-dimer model has a 
simple tensor-network representation on a tilted square lattice. 
As shown in Fig. 1, the partition function TN consists of ten¬ 
sors J 2 ,i 3 ,s 4 located at one half of the plaquettes (the green 
ones) and the diagonal matrices living on the vertices. 
T has four bond indices s,- {i e {1,2,3,4)) corresponding to 
the vertices of the green plaquette. s, e (0,1,2), meaning 
the presence (s, = 1,2) or the absence (s, = 0) of dimers 
on the concerned vertex i,-. In addition, we use i, = 1,2 to 
distinguish different dimers: i, = 1 means a vertical dimer, 
and Si = 2 is a horizontal one. Since the dimers are not 
allowed to touch each other (hard-core condition), there are 
only seven nonzero elements in the plaquette tensor T. The 
nonzero tensor elements T and their corresponding classical 
dimer configurations are schematically shown in Fig. 2 (a-g). 
To, 0 , 0,0 = 1 corresponds to the absence of any dimer on the 
plaquette, Ti^^o.o = To,o,i,i = T 2 , 0 , 2,0 = To, 2 ,o ,2 = 1 represent 
the one-dimer configurations, and 7’i_i_i i(7’2, 2 , 2 , 2 ) = exp(-j6v) 
describes the plaquette with two vertical or horizontal dimers. 
The rest tensor elements are zero (thus forbidden to appear in 
the partition function). 
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FIG. 2: The seven nonzero tensor elements of T and the correspond¬ 
ing dimer configurations on a plaquette. (a) no dimer; (b-e) singlet 
dimer configurations; (f, g) double occupied plaquettes (vertical or 
horizontal). 


The interacting monomer-dimer model under study is de¬ 
fined on a black-white checkerboard lattice schematically de¬ 
picted in Fig. 1 (a), where the dimers are located on the links 
and occupy two lattice sites. Summing over all possible dimer 
coverings, we have the (grand) partition function as 

H = ^ exp [-yS(vVrf - (1) 

k) 

where (c) means the set of all dimer configurations, counts 
the number of the doubly occupied green plaquettes (i.e., there 
are two parallel dimers on each plaquette). v is the coupling 
strength between parallel dimers: v < 0 is an attractive inter¬ 
action, while V > 0 means a repulsive one. In the follows, we 
only consider the attractive case, and set v = -1 as the energy 
scale if not otherwise specified. A,o, is the total dimer num¬ 
ber on the lattice, and ji is the chemical potential of dimers. 


In the fully-packed dimer case, to ensure that every vertex 
is occupied by one (and only one) dimer, a 3 x 3 matrix Ps^x- 
is defined on every vertex, with elements Pq \ - Px q - Po ^2 = 
Pifi = 1. otherwise 0. In the monomer-dimer case, we al¬ 
low monomer doping in the model by setting Fq.o = 1 and 
Po,i - Pi,o - Pq ,2 - Pifl = exp(y6ju/2), with fj. the chemical 
potential. When ^ 00 , the monomer-dimer model recovers 
the fully-packed dimer case. Networking the plaquette ten¬ 
sors T and vertex matrices P, we thus obtain a tensor network 
which faithfully represents the (grand) partition function of 
the interacting monomer-dimer model. 

To calculate the thermodynamics, we adopt the infinite time 
evolution block decimation (iTFBD) method^^’^® for the ac¬ 
curate contraction of the partition function TN. iTFBD was 
proposed for efficient simulations of the time evolution and 
the ground state property (through imaginary-time evolution) 
of ID quantum systems, and then generalized to calculate the 
thermodynamics of 2D classical statistical models^® and ID 
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FIG. 3: (Color online) The specific heat Cy of fully-packed dimer 
model with dimer-dimer interaction v = -1. The results are shown 
to be well converged with various retained bond states 

quantum lattice models.Within the boundary MPS frame¬ 
work, we utilize a kind of “power” method to determine the 
dominating eigenvector (an MPS) of the transfer operator in 
the TN, which consists of a column of T tensors organized 
in a matrix-product operator (MPO) form. However, unlike 
the ordinary power method for matrices, in the transfer MPO 
case the MPS is enlarged after each contraction step, with a 
composite bond space of a direct product of MPS and MPO 
bond bases. Thus the MPS bond dimension grows up expo¬ 
nentially with contraction steps. Therefore, one has to per¬ 
form truncations on the bond space of the enlarged boundary 
MPS, and bring the bond dimension of MPS back to Dc, mak¬ 
ing the contraction procedure sustainable. Performing con¬ 
traction and truncation processes iteratively until the boundary 
MPS converges, we thus obtain the dominating eigenvectors 
of the transfer MPO, with which we can then evaluate the ex¬ 
pectation values of the local observables including the energy, 
the dimer occupation numbers, and the two-point correlation 
functions like dimer-dimer correlations. 

In our practical calculations, we perform the contraction of 
MPS with transfer MPO until the prescribed convergence cri¬ 
terion is reached, say, free energy per site converges to 10“'^ 
(in some cases even down to machine precision). The total 
number of iterations ranges between 4000 and 10^, depending 
on the temperatures and the physical parameters of the model. 
The retained bond dimension of the boundary MPS ~ 100 
~ 150, the convergence with Dc is always checked, the trun¬ 
cation error is less than 10 ® at the critical point, and reaches 
the machine precision (10 away from the critical points. 


III. RESULTS AND DISCUSSIONS 
A. fully-packed dimer model 

Firstly, we investigate the interacting fully-packed dimer 
model on the checkerboard lattice. The specific heat Cy curve 
is shown in Fig. 3, which is computed by taking first-order 



FIG. 4: (Color online) The expectation value (itb) of the dimer 
occupation number on the green plaquettes A (B). The mean value 
n = («x + fig)!! is also plotted. The inset shows the order parameter 
\nA - n\. 


derivative (versus temperature T) of the energy per site. The 
latter is obtained by contracting the TN with one T tensor re¬ 
placed with an impurity tensor T' - v(7’i jq_i + T 2 X 2 , 2 )- From 
Fig. 3, we observe no singularities in the Cy curve, suggest¬ 
ing the absence of any second-order (or lower-order) phase 
transition. 

However, by checking the local occupation number of the 
dimers on the green plaquettes (i.e., dimer density) in Fig. 4, 
we see different dimer densities between the A (dark) 

and B (light) green plaquettes [Fig. 1 (a)] at low temperature 
T <Tc~lA. Especially in the limit T —> 0, the A plaquette 
is filled with a pair of dimers (n^ = 2) while the B plaquette 
is vacant (ng = 0). We use n = (n^ + nB)/2 to denote the 
average dimer number on the green plaquettes. n is verified to 
be a constant in the whole temperature region, this is because 
every site is linked to one dimer in the fully-packed case, and 
each green plaquette contains two sites in net. In the inset of 
Fig. 4, we show that the difference In^ - n\ is nonzero below 
the critical temperature Tc, and it vanishes for T > Tc- There¬ 
fore, the particle number difference Iha “ «l can be regarded as 
an order parameter detecting the phase transition between the 
low-T symmetry breaking phase (n^ 4 n^) and a high-T dis¬ 
ordered phase {nA-ng-l). Because the derivatives of energy 
are always continuous, this phase transition should belong to 
a KT type.^^ 

In order to understand this KT phase transition, we calcu¬ 
lated the correlation length f via the following formula, 

^=l/ln(^), (2) 

where Ai (A 2 ) is the largest (second-largest) eigenvalue of the 
transfer matrix M, obtained by sandwiching two boundary 
MPS tensors 

( 3 ) 

m 

where A is the MPS tensor, with a,b the geometric indices 
and m the sharing physical index. The results are shown in 
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FIG. 5: (Color online) The correlation length ^ versus temperature T. 
The correlations are measured between the green plaquettes of Fig. 1 , 
and the length unit is the distance between two nearest-neighboring 
plaquettes A (or B), i.e., twice the length unit of the original lattice. 
Upper inset shows 1/^ versus 1 / ^jTc — T in the vicinity of Tc (T < 
Tc), with the solid line an exponential fitting. The lower inset shows 
^ versus at T = 0.6 < T,. and T = l.l > T^, the former diverges 
while the latter saturates in the large Dc limit. 



FIG. 6: (Color online) Log-log plot of the correlation function F(L) 
with V = 0 (or, equivalently T = oo limit) for the fully-packed dimer 
case. The fit F(L) ~ reveals the algebraic decaying behavior 
(with Tj obtained by linear fit shown in the inset). 


Fig. 5, where ^ converges rapidly (with increasing Dc) to a 
finite value for T < Tc (lower inset of Fig. 5), verifying the 
existence of a non-critical phase with finite correlation length. 
On the contrary, ^ grows almost linearly with the increase of 
Dc for T > Tc (lower inset of Fig. 5). Therefore, we expect 
the T > Tc region (shaded region in Fig. 5) is a critical phase 
with divergent correlation length. The upper inset of Fig. 5 
is the correlation length ^ as a function of 1 / y/Tc - T when 
T approaches Tc from below {T < Tc). It indicates that the 
correlation length ^ diverges as an exponential of 1 / VIT - Td- 
In Figs. 6 and 7, we present explicitly the dimer number 
correlation function 

r,',; = {OiOj}, (4) 

where = (n^), - n. In Fig. 6, we set v = 0, and the par¬ 
tition function is an equal weight superposition of all possi¬ 


FIG. 7: (Color online) The correlation function at different temper¬ 
atures T (> Tc) for the fully-packed dimer model. Inset: the critical 
exponent rj versus T, rj is extracted from the fit (F ~ L^''). 


ble fully-packed dimer configurations (corresponding to the 
T - oa limit). A log-log plot of F versus L and its algebraic 
fit (F(.r) ~ are shown in Fig. 6, where the algebraic de¬ 
cay is clearly verified. By linear fitting, the exponent rj can be 
obtained, and is shown in Fig. 6. The deviations of F(L) from 
the algebraic behavior (due to numerical errors) can be con¬ 
tinuously corrected by increasing Dc- Notice that in Ref. 13, 
a related exponent is determined to be 77 = 2 for both longitu¬ 
dinal and odd transverse dimer-dimer correlations (and 77 = 4 
for even transverse correlations), from analytical results. Our 
value 77 ^ 1.83 is obtained from dimer occupation number 
correlations Eq. 4 (horizontal and vertical dimers are not dis¬ 
tinguished), which nevertheless well agrees with theirs. 

Besides the T = 00 (y = 0) limit, we also studied corre¬ 
lations of other points in the critical phase with nonzero v, 
and present the results in Fig. 7, from which we can observe 
algebraic decaying behaviors of F for every T > Tc- In the 
inset of Fig. 7, we show that the critical exponent 77 grows 
monotonously with increasing T. 

Therefore, by studying the dimer-dimer correlation func¬ 
tion and correlation length, we find the phase transition oc¬ 
curring at Tc is between an ordered (dimer crystal) phase and 
a critical (algebraic liquid) phase. This scenario also perfectly 
supports the conclusion of a KT-transition at Tc- 

To gain further insight into the underlying physics of this 
critical phase in the fully-packed dimer model, we also ex¬ 
tract the conformal central charge of this system in the critical 
regime T > Tc- The conformal field theory (CFT) tells us 
the conformal invariance at the critical point, and sets useful 
constraints on the critical behaviors of two-dimensional clas¬ 
sical or ID quantum systems^*. The universality class can be 
characterized by the conformal anomaly or central charge c 
of the Virasoro algebra. We use MPS-based method to cal¬ 
culate the central charge, by fitting the block entanglement 
entropy S versus the block size L. Depending on whether 
the system is in a critical or a noncritical regime, the block 
entanglement entropy has different asymptotic behaviors^®’^®. 
In noncritical regimes, S grows monotonously with L before 
saturation; while in critical regimes, the CFT predicts a loga- 
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FIG. 8: (Color online) The block entanglement entropy S as a func¬ 
tion of length L calculated by the iTEBD for the fully-packed dimer 
model, with - 40,60, 80,100. The solid line is the fitting curve to 
Eq. (5) (obtained by linear fit shown in the inset), giving the result of 
central charge c - 1. The inset: the entanglement entropy S versus 
log2(L). 

rithmic divergence^' 

5^^1og2(L) + ^, (5) 

where L measures the site number of the block embedding in 
an infinite MPS, c is the central charge and k is a non-universal 
constant. 

The entanglement entropy S is defined by 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
T 


EIG. 9: (Color online) The specific heat Cy of the monomer-dimer 
model with a chemical potential /j = 0. 



T 


5 =-Tr(plog 2 p), (6) 

where p is the reduced density matrix (DM) of system and can 
be calculated from the converged MPS. However, notice that 
for any L, the dimension of the reduced DM supported by the 
MPS is DI X D^. Therefore it is not possible to capture the 
entanglement entropies for extremely long L; however, by in¬ 
creasing the Dc we are able to simulate the logarithmic diver¬ 
gence for sufficiently long L. By fitting our numerical results 
to Eq. (5), as shown in Fig. 8, we find c = 1 for the fully- 
packed dimer model, in accordance with the result obtained 
by another independent method in Ref. 12,13. 


B. monomer-dimer model 

In this part, we study the interacting monomer-dimer model 
(with p < oo) on the checkerboard lattice. Fig. 9 shows the 
calculated specific heat Cy of the case p - 0, where a diver¬ 
gent peak of Cv occurs at Tc = 0.35, uncovering the existence 
of a second-order phase transition. 

The dimer occupation numbers («b) on the plaquettes A 
(B) and the average n - (n^ + nB)l2 are shown in Fig. 10. 
For T < Tc, the symmetry between the A and B plaque¬ 
ttes is broken { 11 ^ + ng), while for T > Tc this symmetry 
is recovered (ha = vib). In contrast to the fully-packed case 


FIG. 10: (Color online) The dimer occupation number on the plaque¬ 
ttes A(B) and the average value u for p = 0, v = -1. 


in - 1 as a constant), n decreases with increasing temper¬ 
atures in the monomer-dimer case. The limit T ^ co (or, 
equivalently p = v = 0) is an interesting special case, i.e., the 
conventional (noninteracting) monomer-dimer model. The 
mean value n is determined as nr^oo = 0.638 123 109 228 547, 
which agrees perfectly with the previous studies (0.638 12311 
in Ref. 9, and 0.638 1231 in Ref. 10), and provides 15 very 
well converged (correct) digits. The corresponding free en¬ 
ergy per site (negative of the monomer-dimer constant /i 2 ) is 
/ = -0.662798972 833 746 with 15 converged (correct) dig¬ 
its, again in perfect agreements with previous results (/j 2 =- 
0.662 798 972 834 in Ref. 10, and /j 2 =0.662 798 972 7 + 
0.000 000 000 1 in Ref. 11). 

In Fig. 11, we show the correlation length which also 
shows a divergent peak at Tc, the second-order phase transi¬ 
tion point. Notice that in the T > Tc region, the correlation 
length ^ is finite, in contrast to the fully-packed case where f 
diverges in the high-T disordered phase. 

Entropy is another interesting quantity. Actually, we refer 
to two kinds of entropies in the calculations, i.e., the conven¬ 
tional thermodynamic entropy S - iU - F)IT and the for¬ 
mal “entanglement entropy” S e evaluated from the boundary 
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FIG. 11; (Color online) The correlation length for ^ = 0, v = -1. 
The heights of the peaks at grow with the increase of D^. 



T 

FIG. 12: (Color online) The entropy 5 and the entanglement entropy 
5^ for ^ = 0. Because there are A and B plaquettes, the bound¬ 
ary MPS is of period two, leading to even and odd cut entropies S e 
even and odd. Inset: the differential of entropy S for temperature T 
dSIdT. 


MPS. Given the boundary MPS, we can take a Schmidt de¬ 
composition (once for all bonds) of the translation-invariant 
MPS and formally calculate its “entanglement properties”. 
Notice that this bipartite entanglement entropy is different the 
block entanglement entropy discussed above, because the for¬ 
mer is between two half-inhnite chain. As shown in Fig. 12, 
the bipartite entanglement entropy S e shows a clearly diver¬ 
gent peak at Tc, indicating the occurrence of a phase transi¬ 
tion. This observation is quite remarkable, because the con¬ 
ventional thermodynamic entropy S is smooth around Tc, and 
its singularity can only be seen after taking a derivative over T 
(inset of Fig. 12), owing to ||r = Therefore, this “entan¬ 
glement entropy” 5 £ is found to be more sensitive to the phase 
transition (than the thermodynamic entropy S ), and thus can 
serve as an useful numerical tool detecting continuous phase 
transitions. 

In Fig. 13, we show the semi-log plot of the correlation 
functions F for fi - Q. The linear fittings are performed using 
the correlation length f estimated from the transfer matrices 



FIG. 13: (Color online) The correlation function F for different tem¬ 
peratures higher than at p = 0 and their fitting. The fitting used 
correlation length ^ calculated directly (see details in the text). The 
fitting formula is F(L) ~ . 



FIG. 14: (Color online) The phase diagram of the monomer-dimer 
model on the checkerboard lattice with v = -1. The (red) vertical 
line (p = oo) is a critical line; while the blue curved line (with finite 
p) is a second order phase transition line. Inset: Amplification of the 
negative p region, the star denotes the terminating point p^ « -0.51 
of the critical line. 


(Eq. (2)). Note that for both T > Tc and < Tc cases, the cor¬ 
relation functions are exponentially decaying, indicating that 
the high-T phase is non-critical under the monomer doping. 


C. phase diagram 

As a summary of the previous studies of the phase tran¬ 
sitions, we show the phase diagrams of the monomer-dimer 
model in Figs. 14 and 15. The n - T phase diagram (with 
fixed V = -1) is shown in Fig. 14. The red vertical line at 
p = oo is a line consisting of critical points, i.e., a critical line, 
and the KT-transition point Tc ~ 1.4 separates the low-T sym- 
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FIG. 15: (Color online) The phase diagrams of the fully-packed 
dimer and the monomer-dimer (jj = 0) models on the checkerboard 
lattice, (a) The low-T symmetry breaking phase and the high-T criti¬ 
cal phase are separated by a KT phase transition line, (b) The low-T 
symmetry breaking and the high-T disordered phases are separated 
by a second-order phase transition line. 


metry breaking phase and the high-T critical phase. When /j 
is finite, the phase boundary (blue curved) line represents con¬ 
tinuous (second-order) phase transitions, separating the low- 
T ordered and the high-T disordered non-critical phases. The 
blue curved line terminates at jUc, which is denoted by a blue 
star in the inset of Fig. 14. We estimate, by a polynomial 
fitting, that ~ -0.51, below which the low-T symmetry 
breaking phase disappears. 

The v-T phase diagrams of the fully-packed dimer and the 
monomer-dimer (/r = 0) cases on the checkerboard lattice are 
shown in Fig. 15. The phase boundary line of the fully-packed 
dimer model is a KT phase transition line, which separates the 
low-T ordered and the high-T critical phases. On the other 
hand, the phase boundary in the monomer-dimer model with 
ju = 0 is a second-order phase transition line, which sepa¬ 
rates the low-T ordered and the high-T disordered non-critical 
phases. Notice that in both cases, the transition temperature Tc 
vanishes when v = 0, in agreement with the observation that 
the low-T symmetry breaking phase is induced by the dimer- 
dimer attractive interactions v. 


IV. CONCLUSION AND OUTLOOK 

By employing the accurate tensor network method, we have 
systematically studied the interacting monomer-dimer model 
on the checkerboard lattice. The specific heat Cy and the or¬ 
der parameter In^ - n\ show that KT phase transitions occur in 
the interacting fully-packed dimer model (p = oo), in contrast 
to the finite-yu case where second-order phase transitions take 
place. Collecting the phase transition points, we obtain the 
fi-T and v-T phase diagrams with fixed v = -1 andp = 0 or 
oo, respectively. From the phase diagrams, we find that the at¬ 
tractive interactions v < 0 always induce a symmetry breaking 
phase at low temperatures, no matter in the fully-packed case 
ifi = oo) or the monomer-dimer case (fi < oo). Previously, 
people have found similar conclusions for the square-lattice 
interacting dimer models*^. Flere we show that even switching 
off the interactions on one half of the plaquettes (thus reduc¬ 
ing to a checkerboard lattice model), there is still a symmetry 
breaking dimer crystal phase at a low-T. As a consequence, 
the dimer crystal does not break the 90 degree lattice rota¬ 
tional symmetry on the checkerboard lattice. 

The efficient tensor network technique enables us to cal¬ 
culate the thermodynamic properties of the monomer-dimer 
models with a very high precision. For example, the 
monomer-dimer constant can be determined to the machine 
precision. The tensor network method also provides novel 
tools (for example, boundary MPS entanglement entropy) for 
detecting the phase transitions. 

Besides the square and checkerboard lattices, it calls for 
more investigations of this interacting monomer-dimer mod¬ 
els on other lattices, say kagome or star lattice, to explore 
the dimer-dimer interaction effects in more general situations. 
The tensor network method is also applicable for investiga¬ 
tion of these lattice dimer models and we will discuss them 
elsewhere. 
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